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Abstract. A stationary theory of quantum stochastic processes of second or- 
der is outlined. It includes KMS processes in wide sense like the equilibrium 
finite temperature quantum noise given by the Planck's spectral formula. It 
is shown that for each stationary noise there exists a natural output process 
which is identical to the noise in the infinite temperature limit, and flipping 
with the noise if the time is reversed at finite temperature. A canonical Hilbert 
space representation of the quantum noise and the fundamental output process 
is established and a decomposition of their spectra is found. A brief expla- 
nation of quantum stochastic integration with respect to the input-output 
processes is given using only correlation functions. This provides a mathe- 
matical foundation for linear stationary filtering transformations of quantum 
stochastic processes. It is proved that the colored quantum stationary noise 
and its time-reversed version can be obtained in the second order theory by 
a linear nonadapted filtering of the standard vacuum noise uniquely defined 
by the canonical creation and annihilation operators on the spectrum of the 
input-output pair. 



1. Introduction 

In this paper we develop a correlation theory of stationary quantum noise, give 
its spectral analysis and classification and extend the theory of quantum stochastic 
integration PQ to colored quantum noise [2J. The typical example of such noise 
is given by the Planck spectral formula. It shows that finite temperature equilib- 
rium quantum noise is not white in contrast to classical equilibrium noise given by 
the Nyquist spectral formula. The weak coupling limit and rotating wave approx- 
imation [2 0] make it possible to approximate the equilibrium quantum noise by 
canonical pairs of noncommutative white noises in a narrow spectral band near a 
resonant frequency. 

Although such an approximation is sufficient in quantum optics and in many 
other practical cases, it is not satisfactory from a purely theoretical point of view 
because it does not correctly predict the simplest quantum stochastic motion of a 
free Brownian particle in an equilibrium environment, when the resonant frequency 
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In order to describe the output stochastic processes of quantum filters and other 
devices of quantum measurement and communications as functional transforma- 
tions of the equilibrium noise, we also need a generalization of the notion of quan- 
tum output fields |5J |H] in the framework of stochastic integration with respect 
to the colored quantum noise 

The well developed quantum stochastic calculus and noncommutative extension 
of Ito integration to vacuum noise with canonical commutation relations in the time 
domain jS] El is not directly applicable to such cases, neither is the calculus and 
integration with respect to stochastic models of temperature quantum noise with 
flat spectrum 

But we can use the methods of quantum stochastic integration in the spectral 
domain, in which the stationary quantum noise is 5-correlated with a certain mod- 
ification. The latter is required because the colored quantum processes are not 
frequency-stationary and need not be adapted in the spectral representation. Al- 
though non-stationary and non-adapted theories of quantum stochastic calculus 
and integration have been already established |9l ll()llll] . here we will use the much 
more simple approach of mean square integration which is sufficient for the non- 
adapted linear transformations. The corresponding classical theory of stochastic 
integration was developed by Wiener before the Ito integration theory and is based 
on the possibility of representing any (colored) process stationary in the wide sense 
as a linear filtering integration of the standard white noise. 

Here we shall prove that any quantum process stationary in the wide sense can 
be also obtained from a standard one also by filtering. In fact, we will show that 
it can be obtained by mean square integration with respect to a canonical pair of 
orthogonal quantum integrators describing the standard zero temperature (vacuum) 
noise in the second order. Although such a possibility is known for the Gaussian 
case as a Bogolubov transformation which doubles degrees of freedom of the noise 
in such representation, it has not been realized in the theory of quantum noise and 
integration. 

We will remove this unpleasant feature by deriving the fundamental output pro- 
cess for a given quantum noise as the canonically time-reversed noise, which com- 
mutes with the input noise and has maximal possible correlations with it. In the 
classical (or infinite temperature) case, the fundamental output process coincides 
with the noise. In the quantum case it gives the best possible nondemolition filtering 
and time continuous indirect observation of the noise. 

The fundamental output process for a given Gaussian quantum noise in the 
framework of quantum stochastic calculus was introduced in |12| and the possi- 
bility of its nondemolition observation was demonstrated for Markovian models of 
interaction with a quantum open system in the quantum theory of filtering (see the 
recent survey and papers cited therein). 

To explain the idea of the fundamental output process, let us represent a canon- 
ical Bose-noise annihilation and creation pair b, b^ with non-zero temperature cor- 
relations 

(b f b)=n, (bb t )=n+l 

by the linear combination b = \Jn + la + yfnc of the zero temperature pairs a, a^ 
and c, c^ with 
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[a,at] = l, [c+,c] = l, [a f ,c] = [a,c] = 0. 

Then the combination b = y/na, + \fn + lc defines a canonical pair b, b^ of output 
creation and annihilation operators with correlations 

(b t b)=n+l, (bb t )=n 

which commutes with the noise, [b,b] = [b, b 1 ] = 0, and provides the maximal 

mutual correlations (bb^) = ^/n(n + 1), corresponding to correlation coefficient 
r=l. 

In the classical case, the latter would mean that b = b, but this is not the case: 
b commutes with b^ but b does not. Nevertheless b is uniquely defined for the 
noise b by the described properties as the time-reversed noise, which in the real 
arrow of time plays the role of the fundamental output process. Moreover, the zero 
temperature pair a, c is uniquely determined by the pair (b, b), being given by 

a = vn + f b — y^nb, c = Vn+I b — \/nh. 

2. Quantum Correlations and Reversed Processes 

Here we sketch of a second order theory of quantum noise and linear filtering 
of stationary signals, initiated more than 20 years ago in the pioneer paper 14 . 
A real scalar quantum stochastic noise as a process x. = (xy) in second order is 
completely determined by the zero mean values (xj) = and finite variances alone 
(xj) < oo of the Hermitian variables Xj = xj and also by their not necessarily 
symmetric covariances, or correlations (xjXj) ^ (xjXi). The products XiXj in (•) 
form a Hermitian-positive kernels (xi\xj) which define the scalar products in a 
complex Euclidean space. By writing an invertible correlation kernel (x^Xj) = 
(xj|xj) as a bra-ket inner product (xj| ■ \xj), we can describe the quantum noise in 
second order by real vectors \xj) — |xj)" equipped with the involution |z) i— » |z)" = 
|zt). The latter, defined by 

lEO x i)*=ECW VC.eC 

such that |z)"" = |z), gives vector representation of the Hermitian conjugation 
z ^ z * = E C j x j with Cj = Cj f° r the complex combinations z = Cj x j = C ' x °f 

the Hermitian variables x^ given by raws C. = (Ci)- 

Let us interpret the index j as "discrete time", e.g. j £ Z when x. = (xj) jgZ 
is two sided infinite sequence (■ • • , x_i, xq, xi, ■ • • ) of Hermitian variables with the 
time ordered correlations 

Kij = (xi xj) = (xi\xj) , i > j e Z 

and the reverse ordered Kji = (xj\xi) which we need to define the scalar product of 
the corresponding vectors |xj). In contrary to the classical theory, these correlations 
need not be real valued even for Hermitian Xj , so that time reversal corresponds to 
their complex conjugation: 



Kji — (xj Xi) — ( Xi Xj) — Kij. 
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Although the reversed correlations may be considered as not "observable" in the 
real arrow of time, they are defined by complex conjugation in such a way that the 
Hermitian matrix K = [Kij] is positive definite: 

as is the complex conjugate = transpose K — [Kij] — [Kjj] = K : 

We claim that K can be treated as the correlation matrix 

= (xj Xj) = (xi|x 3 ), Vi > j £ Z 

of another sequence Xj = (• • • , x_i,xo,xi, • • • ) of Hermitian x, commuting with Xj 
for all j £ Z which is maximally correlated to x. = (xj) je z m the sense that the real 
covariance matrix G = [(xj Xj)] ispositive and symmetric as the geometric mean 

/ ~\V2 _ 
G = K ■ K) of K and K = K. The sequence x describes the fundamental 

output process as reversed noise by the quantum Hermitian variables Xj = xj 
which will also be represented by real vectors \x.j) = |xj) b with respect to another 
involution |z) i— > Iz) 11 uniquely defined on the complex combination =C ■ x by 
complex conjugation of the coefficients ( ; e C. This suggests a natural norm 

||C.|| 2 = (zt z ) + (z J) = c(k + k) -c 

in the space E of test sequences £. = (Q) for the simultaneous treatment of linear 
combinations \z) = ^2Cj\ x j) °f t ne i _ rea l vectors |xj) = |xj)", and the combina- 
tions |z) = Y^Cjfx-j) of b - real vectors |x,-) = |xj) b , representing Xj = xj. 

From now on for the sake of simplicity we shall restrict ourselves to the case 
when the matrices K and K commute as it is in the stationary case 

i^'ij — ki-j ■ k — j - — - kj , ^3 ^ ^ 

(see also ^21 for a more general, noncommutative case.) In this case that covari- 
ances of x.; and Xj are described by the real symmetric matrix 

G=(KK) 1/2 = (K K) 1/2 = G 

as geometric mean of the commuting K and corresponding to maximal correla- 
tions with zero commutators between x and x: 

Gij = (xj Xj) = (xjXi) = Gj h Vi > j £ Z 

If K is invertible, this output process is identically correlated with the noise, corre- 
sponding to the correlation coefficient r — I , and in the classical case K — K such 
the sequence (Xj) is always identified with (xj). The invertible case corresponds to 
the thermal noise which we call standard if Gij — Sij , that is if the inverse noise 
fx.,) is (^-correlated with (xj) such that K = K^ 1 . The opposite situation KK = 
when (xj) is uncorrelated with (xj) such that r — is possible only in the quantum 
case. It corresponds to the vacuum noise which we call also standard if K + K = I. 
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Let us prove the existence of the time-reversed sequence x / x in the second 
order quantum theory corresponding to the case K ^ K, and find the adjoint 
involution \z) i— ► |z) b = |zt), 

l£C^> = £C*l^) VC,eC 

such that |z) bb = |z), giving the reversed vector representation of the Hermitian 
conjugation z i— »zt = ^ £*Xj for the complex combinations z = ^ ^Xj in the case 
of invertible X. 

The commuting matrices K and if have common eigen vectors u. (generalized 
row- vectors) , given as complex eigen-sequences u. = {uj)- eZ such that 

E = ' X/ = kuj, j E Z 

ieZ ie z 

with eigenvalues k and k > 0, corresponding to the complex conjugate eigen se- 
quences TT. = (uj)j €Z : 

*^^HiKij = KUj, ^^viiKij = kuj, j e Z 
iez iez 

Indeed, let K and K have a common eigen-row- vector u. with eigenvalues k > 
and k > respectively. Because if is the complex conjugate matrix of K, TT. is 
also a common eigenvector, but with the real eigenvalue k for K and k for K. Two 
cases arise 

k = R: In this case we may assume without loss of generality that Uj = 
uj, Vj g Z 

k: In this case the eigenvectors u., u~. are orthogonal and so we have 

This means that the index set SI of the eigen sequences (v\ := {uj(v))j e z is 
equipped with a natural flip 

v i ► — i/, —{—v) = v, = k(—v) 

such that u~.(v) = u. {—v) and the (generalized) Plansherel measure dv on SI is invari- 
ant under this flip. The complete (with respect to dv) orthogonal set {u.(v)\v e SI} 

can be chosen in such a way that Uj(v) = Uj(—v) iff k(u) = k(—u) and ^ u?(v) = 

iez 

iff k(v) 7^ k(— f). This flip defined an isometric involution p* (y) := Tp {—v) = tp* (v) 
in the Hilbert space L 2 (SI) as antilincar map p i— > p*, p** = p satisfying 

(</>*|V) = J V {-v) tp{v)dv = j p (i/) ^ (-i/) Av = (i/>*\<p) \/<p, ip e £ 2 («) . 

We shall assume without loss of generality that the Plansherel measure dv on 
the null subset N = {V : k(z/) — = k(^)} is zero, taking N = such that its 
support SI is identified with the union U N^: , where 

= {v : k{v) > 0} , N± = {v : k{v) > 0} . 

are the complementary subsets of the null sets N + and N_ for k e L 1 (N+) 
and k e L 1 (N^). Let us prove that the pairs {SI 3 v Xj(v)} G L 2 (NiJ:), 
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{!)9i/h Xj(u)} G L 2 (Nf) of complex amplitudes 

Xj{u) = k{v) 1/2 u 3 {v) = (i/| X j), Xj(u) = k{v) 1/2 u 3 (v) = 
which are related by the isometric involution * as 

Xj = ^{—v^^Ujiy) = ^{—v^^u^—v) = Xj 

where u* = u, give the spectral representation of the quantum noise vectors \xj) 
and |xj). 

Indeed, by virtue of the completeness of {u.{v) \ v G 0} it follows that 
x\xj := / Xi(u)*Xj(i/)dv = / Ti{v)ui{v)uj{v)ilv — Kij 



x\dtj :— j Xi(v)* Xj{v)Av = j ~f{v)ui{v)uj{v)&v = dj 

where we used the invariance of the measure on f2 with respect to the flip v — v. 
Let us define the antilinear maps |z) i— > |z)", |z) i— » |z) b of vector conjugation 

z^A^/V, z^A^z*, 

respectively on the dense domains 2?j C L 2 (N|:) and 2?|, C L 2 (Nf ) as the subspaces 

P B = {z: A 1/2 z G L 2 (^)} ,X> b = {z : A~ 1/2 z e £ 2 (f!)} 

of square integrable complex amplitudes z G L 2 (0) with the support in N;j: and 
z E L 2 (O) with the support in where A 1 / 2 = y/WJ k and A -1 / 2 = J k/k are 

~ 1 /2 

well defined as positive operators of multiplications respectively by [k(v)/k(u)] ' 
and [k{v)/k(v)\ 1 / 2 into the subspace L 2 (6) the common support 

Q = {v:k{v)^{)^ k(v)} = N| n Nf. 

Note that 2?j C L 2 (f2) is generated by the complex amplitudes z = ^2 (jXj with 
z* G T>\, and the adjoint conjugation on the subspace T>\, C L 2 (f2) is generated by 
z = ECj^j with z* G 2?(t such that ^tt = ^b- Moreover, z^* = i* b , z b * = z*» and 
z b = Az", z" = A _1 z b are densely defined involutions in the subspace L 2 (6) and 
A* = A" 1 on 6. 

In the case of nonzero temperature, when 9 = and Dj, 2?t> are dense in L 2 (O), 
the spectral representations of quantum noise and the output process are connected 
by the spectral linear filters 



and X, — Xj , x^ — Xj for the generating spectral sequences (xj ) and (£j ) : 



x{v) = A(z^) 1 ^ 2 'x{y), x{v) = X(v) 1 ^ 2 x(i>), 

x){u) = \{vy 1/2 Xj {v) = x 3 {u), v G N±, 
x)(v) = \{v) l/2 x 3 {v) = Xj{y), i/eNj. 
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The latter implies |zt) = |z)", |z^) = |z) b for the complex-linear combinations z 
ECjXi, z = £Ci5q since 
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The spectral amplitudes the time representation 

Xj 1 = j Ui{v)xj(v)dv , Xj 1 — J Ui(v)xj(v)dv 

describe the canonical vector realization Xj,Xj of \xj), \xj). It is given by the matrix 
elements of the square roots X = K 1 / 2 and X = K x l 2 as 

Xj — , X j — Xi j , G Z 



.. . V 



such that x 1 , = Xj 1 as X = X. They are thought as the £ 2 (Z) columns 



Xj = (x/) , x, - (V) 

with the scalar products x\xj = (xj|xj), x\xj — (xj|xj): 

X\} Xj ^ ^ X kiXfcj -^ij ^ ^ XfciX fcj Xj X^ ; X\}xj G 'ij . 

fcez fcez 
Hence the inversion x 1— > x is represented by the usual complex conjugation x = 
x, Xj = ( Xj l ) ie , which coincides with transposition : Xj 1 = xj due to the 

sclfadjointncss of the square roots of K = and so K = K . Although the 
vectors Xj 7 Xj seem to be complex if Xj ^ Xj due to Xj := xj = Xj, they are 
selfadjoint x j ^ — x j 5 x j — x j in the case = f2 with respect to the involution 

-J = L-^z*, z° = L^ 2 z*, where z* = z: 

1/ X j — ^ Ij feXj - — - Xj ^ Ij X j — ^ Ij f^Xj — X j 

are the linear input-output and reversed filters in the time representation. The 
operator L = KK^ 1 with the matrix elements 



Lij — 



j Ui(^)A(^)iij(^)dj/, = k(V)/k(V) 



is characterized by the properties L > 0, L — L 1 and is called the modular operator 
in this invertible case. It can be described by the adjoint involutions as L\z) = |z)'\ 
and its algebraic analogue is the main object of study in Tomita-Takcsaki theory. 
Note that if K is not invertible, the canonical realizations of x and x are still 
connected by complex conjugation as the isometric involution z 1— >z*, although it 
may not be described by the linear filtering because the adjoint involutions (j, b 
may not be represented in £. In contrast to the classical theory the constructed 
canonical realization x 3 of the reversed quantum process does not coincide with Xj 
and is even orthogonal to the real quantum noise if N + U N_ = ft. 
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In the case of a stationary time sequence xj = y/ex(tj), tj — sj the spectrum Q 
is in the interval [— e _1 /2, £~ 1 /2] C R, the flip v — v is the usual reflection, and 

u_j(v) = e 1 / 2 cxp{27ri^^} — uj(iy). 

The canonical vector realizations of the stationary sequences (xj) and ( Xj) is given 
in the time representation by the matrix elements Xj 1 = Xj-i, Xj = Xj-i-> where 
the complex sequences 



Xj = e 1 ^ 2 / k(^) 1//2 cxp{— 2-K\vtj\Av = / x 3 - (y) Av 
Jn Jn 

Xj = e 1 ^ 2 / k(v) 1 ^ 2 cxp{— 2-K\vtj\Av = / Xj(v)Av 
Jn Jn 

are connected by the usual time reflection x(tj) = x(—tj) with respect to t = 0: 
Xj = X-j- If K is invertible, this time reflection is described by the stationary 
input-output and reversed linear complex filters 

X-i — ' ^ i-jXji Xi — I ^ i-jX-ji 

jez jez 
where l x / 2 i-j = L 1 / 2 ^, l~ 1 / 2 i -j = L~ 1//2 y are given by the integrals 

Iq 1 ^ 2 j :— e / A(^) 1 ^ 2 exp{27ri^j}d^, Iq" 1 ! 2 j := e / A(^) -1 ^ 2 exp{27ri^tj}df 
Je Je 

over the support 6 C of 7 = (kk) 1 ^ 2 in the general case. 

3. Spectral Decomposition and Linear Filtering 

Let us consider the stationary quantum noise x:fn x(t) with continuous time 
t G R. It can be treated as the limit as e — > of a stationary sequence e~ 1 ^ 2 x e j 
with correlations k £ j = ek(t), given by a complex positive-definite function 

/oo 
exp{ 2irivt} K(v)dv = (x(t)|x(0)) . 
-00 

This may exist only as a generalized function(distribution) if the support Nzj: C R of 
the spectral density n is unbounded. The reversed noise x : 1 1— > x(i), corresponding 
to the spectral density k(^) = k(— u), is described up to second order by the 
autocorrelation function 



k{t)= I exp{27rii4}K(i/)di/ = (x(i)|x~(0)) 

J —OO 

and by the symmetric cross-correlation function 

exp{ 2-K\vt) 7(i/)dzv = (x(i)|x(0)) , 

-OO 

where 7(1/) = , so that r is the convolutional square root of 

k*k (t)= k(s~ t)k(s)ds. 

J —OO 
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If k * k = 0, i.e. if kk = 0, and so N^: n = wc have stationary vacuum noise. 
The stationary noise is called standard vacuum noise if k(t) + k(t) is the Dirac 
5- function, i.e. if k + k = 1. 

Such a noise is purely nonclassical because the condition kk = for a symmetric 
function k = k is only possible in the trivial case N;j: = = when x = = x. 

In the general case the subset 6 = N^nN^: of il is not empty but the complement 
Q 1 - = £1 \ can be decomposed into the disjoint union of 

N_ = {v e Q\k{v) = 0} C Nj, N+ = {v e CI\k(v) = 0} C Nf, N+ n N_ = 0. 

This means that the quantum noise and its time-reversed version can be uniquely 
decomposed in the correlation theory into sums x c + xe , x D + xe of uncorrelated 
vacuum (x ,x ) and thermal (xe,xe) components. In the spectral representation 
they are given by 

x Q {t) = P±*(t), x (t) = P±x(t), 
x e (t) = Pex(t), x e (t) = P®x{t), Vt e R 
where P+ and are the projectors on L 2 (N + ) and L 2 (N_) and P® = P+P- is 
the orthoprojector on the subspace L 2 (Q). The orthoprojectors P + = 1 — 1 + and 
P_ = 1 — 1_ define the best input-output and output-input linear estimates 

P+x(t) = l®x{t) = \ 1/2 x(t) of x (t) = X @ 1/2 x(t) + x (t) , 

P-x(t) = l&x(t) = \ ~ 1/2 x(t) of x(t) = \ @ - 1/2 x(t) + x (t) 

where le, defined by \&{v) = 1 if f € 0, !®{v) — if v £ 0, is the characteristic 
function of = n N^, 1+ and 1_ are the characteristic functions of N + and 
N_ , and 

AeM 1/2 = (?MM,))l e M. 
In the time representation 

x{t) = {x(t ~ *)|* G R} , x{t) = {x(t - s)\s G R} 

x (t-s)= / K(^) 1/2 exp{27ri^(.s-t)}d^ = x(s-t) 
Jn 

xq (t) and 5?e (i) are obtained by the replacing the interval Q of integration in \ = Xn 
and x = Xn by 0, with x a (t) given by x N _ and x (t) given by x N+ . The optimal 
input-output and output-input filters are given in the Fourier representation by the 
complex stationary linear nonadapted integrals 

/CO 
le 1,2 {s-r)x{t-r)dr^[ll( 2 *x) (s) 
-CO 

/OO 
l e - 1/2 (s - r)x(t - r)dr = (l e 1/2 * x) (s) . 
-oo 

Here Z e 1/2 (s - t) = L e 1/2 (s,t), l e ~ 1/2 (s - t) = L e ^ 1/2 {s,t) are the complex 
positive definite generalized functions 

l e 1/2 (t)= J AeM 1/2 exp{27rn4}d^ l e ~ 1/2 {t) = j \ @ (u)~ 1/2 cxp {2-nwt} Av 

characterized by the modular property 

/e 1/2 (-t)=7e 1/2 (*) = /e- 1/2 W. 
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The case N+ = = N_ corresponds to the purely thermal noise which is called 
standard thermal noise, if kk = 1 in a given spectral region ft C R. If £1 = R, 



it can be written as 



k * k 



(t) = S(t) in terms of the Dirac correlation function 
r(t)=6(t). 

The thermal noise is called white if the spectrum k and thus also k is flat: 
k(v) = a 2 = k(v), W € White noise coincides with its time-reversed version 
and is essentially classical (at least in second order). There exists just one such 
standard noise x = w = x described in the second order theory by the correlation 
function k(t) = S(t) = k(t). 

In the quantum case there are many standard thermal noises and ^-correlated 
with them reversed version \{t) = x(— t). They are parametrized by a modular 
spectral function 

\{v) > o, \{v) = \(-v) = A(^)-\ Vv g n 

giving the standard correlation functions k, k as the convolution square roots 

k(t) = r 1/2 {t), k(t) = i 1/2 {t), vt g r 

Any stationary quantum process y : t — > y(t) and its reverse y : t — > y(t) described 
up to second order by the correlation function 

(y(t)y(O)) = J exp{2^i^}^(^) 2 d^ - (y(O)y(t)) 

can be obtained by stationary filtering of standard quantum noises x and x, i.e. 
such that their spectral densities 

k{v) = 1_{v) + \ @ {v)- 1/2 ,k{v) = + \ B (v) 1/2 

and 7(V) = le(^). In the spectral representation 

y(t) = {y{t, v)\v G ft} , y(t) = {y(t, u)\v G ft} 

the filtering is given by 

i){t, v) = f{v)x{t, v), y(t, v) = f{u)x(t, v) 

where / = / = / is a real symmetric transmission function / : Q — > R + of the 
standard complex amplitudes 

x(t, v) = cxp{27rij/t}x(^), x(t, v) = cxp{27ri^t}x(^) 

x(y) = \_{y) + A e (^)" 1/4 , x{y) = + A e (z/) 1/4 . 

The transmission junction f is uniquely defined in the space L 2 (f2) by 

/(i/) = W.)?(.)] 1/2 leM + [a(iz) V al^)] 1 / 2 !^^). 

In the time representation y(t) = {ip(t — s)}, y(t) = {ip{t — s)}, given by the Fourier 
integral 



ip(t) = J y{v) cxp{27ri^t}d^, ip(t) = J y(f) exp{27ri^t}d^, 
the filters are written as real nonadapted stationary transformations 

/GO ^ 
(j>(s — r)x(t — r)dr = ip(s — t) 
-oo 
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where 

<f>(t) = J f{v) exp{27rii^}di/. 

An important example of the stationary quantum thermal noise is given by the 
Planck's spectral density 

k(v) = hv{^ hv - l)" 1 , k{v) = hv{l - e-P hv )~\ 

In this case k(v) — k(v) = hv > if v > 0, i.e. 6+ = R + , \(y) = cxp{/3/iz/} ^ 
Vv £ R, and A(V) ^ 1 for v ^ 0, and \(v) — 1 as for the classical white noise only 
if p = 0. It describes stationary Bose noise in which the electromagnetic (optical) 
field is in an equilibrium state of the temperature 9 — 1/(3. The modular function 
for such noise, 

l(t) = J exp{2mvt} exp{f3hv}dv = S (t + = CX P jy/ 3 ^ j 
is extremely singular, as is the correlation function 

(h = h/2n) with complex conjugate 

hf. f h„d IV 1 d 



fc (,) = ^l-exp|- T/3 -|j -S { t) 

for the time-reversed Planck noise corresponding to the density k. 
Although the cross correlation function 



rit) = j 



Sm 1 2^ dt 



is less singular, due to the integrability of its spectral density 



7(1/) = fti/[exp{/3/ii//2} - exp{-/?hi//2}] \ 

quantum stochastic integrals corresponding to linear stationary filters for such noise 
are not defined in the usual sense. Indeed, the standard quantum pair (x,x) corre- 
sponding to the Planck density have the singular correlation functions 

l~ l ' 2 it) = Sit + ihp/2), l x ' 2 it) = Sit - ihp/2) 

and the cross-correlation function r(t) = S(t). To obtain the Planck noise y and its 
time-reversed version y one should integrate the standard quantum noise x and x 
with the non singular filtering function 

f(t) = J [^/2sh(/3^/2)] 1/2 exp{27ri^}d^ 

But the result of the integration remains singular, corresponding to the singularity 
of k and k. 
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4. Quantum Stochastic Integration 

In quantum theory the noise x (i) is thought of as represented by Hermitian, or 
even selfadjoint operators in a Hilbert space of state vectors (probability ampli- 
tudes), with zero mean values c/)[x(i)] = (c/)\x(t)\(j>) — 0. Here \cf>) is a fixed unit 
vector of the given state </> which can be always be chosen as vector state represented 
by the vacuum vector | (j)) = 5$ in a Fock space TL = T . The quantum correlations 
K(s,t) = k(s — t) are defined for a stationary state <fi by the scalar product in TL, 

k(t) = (cj> |x(t)x(O)|0) = 0[x(t)x(O)]. 

If k(t) is complex, the operator x (t) =x(i)^ will not commute with x(0) =x(0)^. The 
time-reversed correlations k(t) = k(—t) are "observable" in the real arrow of time 
only for the time-reversed noise x which gives a natural representation in the same 
Hilbert space TL an output process x : t i— > x(t), for which k(t) — (<f> |x(i)x(0)| <j)) 
with the same <f> G TL. 

The vector processes t i— > |x(i)), t i— > |x(i)) are obtained from the operator 
representations in TL by the identifications 

| x(t)) := x(t)| 0) = x (t) , |x (t)) := x(t)| 0) = 2 (t) 

in the minimal subspace generated by the set {x(t)| 0) , x(t)| cj))} C TL. In a classical 
theory this gives a one-to-one correspondence between vectors and operators of 
multiplication in the total space TL. This is not the case in quantum probability 
theory as one can observe for vacuum correlations, for which it is possible that 
|z) = z\4>) =0 for annihilation operators z =/= 0. This is why the vector integration 
theory is not sufficient in the quantum case; we need an operator form of quantum 
stochastic integration. 

According to the nondemolition (causality) principle of quantum theory |12| the 
input observables x(0) must commute with the fundamental output observables 
x(t) for all t < and they must commute also for t > because the noise x can 
be regarded as the output process for the time-reversed noise x. At the correlation 
level this is expressed in terms of the reflection symmetry 

(0 |x(0)x(t)| 0) = r(t) = (0 |x(i)x(0)| 4>) , VieR 

~ 1/2 

for the cross correlation function r(t) = (k * k) (t). 

Now we are going to develop a second order theory of quantum stochastic in- 
tegration with respect to any colored quantum noise and its fundamental output 
process. The realization of such theory was given in jl()j in terms of a representa- 
tion of the canonical commutation relations in the symmetric Fock space TL. This 
provides a one-to-one correspondence between second order integration theory and 
the Gaussian integration theory of quantum thermal noise with respect to the vac- 
uum state 4>. Without loss of generality we may assume that the pare x (t) ,x(t) is 
standard, i.e. kk = 1. 

Let us denote by £ the complex test functions 1 1— > ((t), for which the "integrals" 



/oo />oc 
c(i)x(t)dt, c-x= / cm(t)dt 
-oo J — oo 
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are expected to be well defined as selfadjoint operators in H if C = C If C is given 
as a Fourier integral 



i(t) := J c{v) exp {2'Kivt} dv = c(— t), 



then the operator- valued functionals £(—>£' x iC^C' x can be treated as quantum 
stochastic integrals 

c-x = y c(i/)X(di/), y c(i/)X(di/) = c • x 

for the Fourier transforms c = C of ( = c and c = Q oi C, = c. 

The quantum stochastic integrators X, X, defined in |14| as operator-valued 
measures on the spectral intervals A C £1, in second order are described by spectral 
vector-measures 

|X(A)) := X(A)|0) ^X (A) , |X(A)) := X(A)|0) ^X (A) 

which are related by isometric involution * = — * on the vector space as X* = X, 
i.e. X* (-A) =X(A). 

They, together with X (A) f = (X(A)|, X (A) f = (X (A) | are assumed to satisfy 
the following properties: 

(i) orthogonal cr-additivity for disjoint unions ^2 Ai = A 

X{A) = ®X(Ai), X(A) = ®X(Ai). 

(ii) absolute continuity: 

A > (A) t X(A') = f n{v)dv 
J AnA' 

A > (A) t A > (A') = f k(v)&v 
J AnA' 

A > (A) t A > (A') = f 7 (j/)di/ 
J AnA' 

(iii) selfadjointncss on = N:J: D Nf: 

X(A) t) = A -1 ■ X(A) = X(A), X(Af = A ■ X(A) = X{A) A C 9, 

where 7 = («) 1/2 , X(u) = {k/ K ) 1/2 and A • X (A) = J®\(v)X(dv). The cr- 
additivity makes it possible to define this integral as integral as an orthogonal 
vector measure X absolutely continuous with X on 0. The condition (ii) can be 
symbolically written in the form of a multiplication table 

X(d^X(dv') = 5 vv ,K{ V )d V , X(dv?X(dis') 1 = 5 vvn {v)dv 

X{dv)^X{dv') = 5 vvll {v)dv, X{d^X{dv') = 8 vu .K{u)dv 

where S vv ' = if v 7^ v ,S VI/ / — 1 if v = v , and X (dv) — \{y) X {dv). An 
operator realization of this table can be given in the Fock space T with respect to 
the vacuum state vector | <f>) = Ti, JOj • 

Using this multiplication table we obtain the isometry property 



|yVl 



a(v)K(h') 1 ^ 2 + c(v)n(v) i ^ 2 



2 



'dv = tfb 
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of the map y i— > 6 — an 1 / 2 +ck 1 ^ 2 for the operators y = a-x+c-x = rj-x with rj = d+c 
on the state vector (f> into the spectral representation b(y) = (u | y) of | y) = | r\ ■ x). 
This extends the quantum stochastic integral r\ ■ x with fotfo = J |6(i/)| 2 dz/ < oo 
from the functions 77 (i) with simple (a, c) to any functions n = a + c for which the 
complex amplitude b remains square integrable. To ensure also (<j> |yy'| </>) < 00, 
we must add the condition of square integrability b* S L 2 (il) for the amplitude 
b* = a* k 1 / 2 + c^k 1 / 2 , that is, we define the quantum stochastic integral y together 
with yT only for the test functions a,c C 8. In the classical when 
b* = b*, the square integrability of b and b* are equivalent to that of a and c with 
respect to the symmetric measure k &v — k Av. In the purely quantum case of 
vacuum noise it is not so, and the conditions 6,6* € L 2 (VL) are equivalent to the 
square integrability of a and c with respect to n V ac dzA 

In spectral representation the standard vacuum noise x (t) together with the 
reversed process x (t) is described by the canonical operator- valued measures 

A+ (A) - { ?( A )' A C N + A (A) — { * (A) ' A C N + 

A (Aj "\ X(A),AcNi ' A - {/X >~ { X (A),AcN^ 

of independent creation A A = A + (A) and annihilation A A = A_ (A) such that 

<0|A+=O, A A |0)=O. 

It defines the standard operator vacuum measures X and X as 

X(A) = A_ (An Nf) + A+ (A n N|) , X(A) =A_ (An N j) + A+ (A n N± ) 

on the intervals A C SI for the given state | <j>) G 7i. The canonical pair (A^, A^j~) of 
the creation and annihilation measures is characterized by the following properties 
of fiip-adjointness 

A A t — A+ = A A A +t — A~ A = A~ 

and orthogonality of all products apart of A_(A)A+(A') for An A 7^ 0: 

(0 |A_(-A)A+(A')| 0) = ^(A n A'), (0 |A+(-A)A_(A')| 0) = 0, 

(0 |A+(-A)A+(A')| 0) = 0, (0 |A_(-A)A_(A')| 0) = 0. 
This table can be written also symbolically in the canonical form jS] as 

A-(di/)A+(di/) = 6 uv >dv, A + (di/)A _ (di/) = 0, 

A + (df)A + (dz/) = 0, A_(di/)A~(dz/) = 0, 

for all f, 1/ € SI. This table includes also A~(dz/)A + (dz/) = 5 vlJ 'dv and all other 
products equal zero by the reflection v t— > — v. 

Let us prove that an arbitrary (not necessary vacuum) stationary (in the second 
order sense) quantum process y (t) together with its time-reversed version y (t) can 
be obtained by quantum stochastic integration with respect to the flip-selfadjoint 
canonical pair (A + ,A_) over SI C R, or, equivalently, with respect to the self- 
adjoint quadruple (A + , A~; A_, A+) on the positive part Sl + C R + of SI. 

Indeed, in general the pair 

y(t) = / exp{-2?ri^}Y(d^), y(t) = / exp{-27rii/i} Y(dv) 
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of operator-valued distributions c i— > c • y, c • y is given by a pair (Y, Y) of operator- 
valued orthogonal flip-selfadjoint measures Y (dv)^ = Y (— d^) with multiplication 
table 

Y(di/)Y(di/) f = a(vfdv, Y(di/)Y(d^) f = a(v)a(v)dv, 

Y(di/)Y(dt/) f = a(v)a(v)dv, Y(Av)Y(dv) ] = a(vfdv. 
Such measures can be obtained as quantum stochastic integrals 



Y(A) = / a(v)k-(dv) + a(v)A + (dv) 

J A 



Y(A)= / a(v)A-(dv) + a(v)A+(dv) 



A 



from the canonical pair (A + , A_) of flip-adjoint annihilation and creation measures. 
Moreover, in the nonclassical case a(v) ^ a(v),Vv e ft, corresponding to quantum 
noise of zero or finite temperature, the canonical pair A = (A + ,A_) is uniquely 
defined over f2 by the pair (Y, Y) from 

(a(vf - cr(z/) 2 ) A+(dv) = a(v)Y(dv) - a(v)Y(dv), 

(a(vf - cr(» 2 ) A_(d^) = a(v)Y(dv) - a(v)Y(dv). 

This means that the direct and reversed quantum stochastic integrals £ ■ y, ( ■ y can 
be written for the Fourier integrals of ( = c as 

C • y = / 9 iy) A- (dv) + g (v) A+ (dv) = A(ac, ac) 

with g (v) = a(v)c(v) and 

C • y = J 9 iy) A_ (dv) + g (v) A+ (dv) = A( ac, ac) 

with g (v) = a(v)c(v). 

To obtain these integrals in the time representation as linear stationary filters 
of the adjoint pair a = (a + ,a_) of the canonical annihilation a_ (t) and creation 
a+ (t) distributions 

a_ (t) = J exp {-2mvt} A_ (dv) = a+ (t) 1 , 

a + (i) = y cxp{-27rii4} A+(dv) = a_ , 

or their time reversals a - (t) — a_ (— t), a + (i) = a+ (—t) we can use the Fourier- 
Parseval identity 

A(53~,£+) = a,(tp~,<p + ). 
The latter is given by the standard creation and annihilation integrators 

A+ (dt ) = a+ (t ) dt = A_ (dt) f , A_ (dt ) = a_ (t ) di = A+ (dt ) * 

in the time representation as the standard quantum stochastic integral 



&(ip- 



/OO 
v + (t)k+(dt) + <p- (t)A_(dt) 
-oo 
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Applying this to the quantum stochastic integral in the spectral domain with the 
general integrands 

/" (v) = a(v)a(v) + c(v)a(u) = ^ (v) 

f+ W) = a{v)a{v) + c{v)a{v) = p + (y) 
corresponding to a ■ y + c- y = A(/~, / + ), we obtain a -y + c- y = a( ip~ , <p + ), where 
the components 

<p- ( t ) = (a * + (* * v)(t) = r (*) 

<P+ (*) = * €)(*) + * »/)(-*) = /+ (*) 
are given by the convolutions 



(JO 

(<r * C) (t) = J a(t-r)( (r) dr 



with transition function 

a{t) = J exp {27ri^i} a(v)diy = a (— t) 

Thus, the generalized processes y and y are obtained by the stationary filters 

/oo />oo 
a(s - t)A_(ds) + / (7(s-i)A + (ds) 
-oo J —oo 

/OO /" oc 

CT(s-t)A_(ds)+ / (t(s - i)A+(ds) 
-OO J —oo 

corresponding to the quantum stochastic integral representation 

/CO pOO ^ 

(a*0(s)A_(ds)+ / (<7*C)( S )A+(d S ), 
-oo J —oo 

/OO pOO 
(a*()( S )A_(ds)+ (a*C)( S )A+(d,s) 
-oo J — oo 

for the test functions ((t) = ((—t) with the square integrable (a and (a. 

Conclusion 

Thus we proved that for each quantum stationary noise there exists a fundamen- 
tal output process which replaces the noise if the time is reversed. The nondemo- 
lition observation of the quantum noise via the measurement of the fundamental 
output process provides the best mean square filtering of the noise and in the clas- 
sical limit completely eliminates this noise. 

The input and output stationary processes can be decomposed in the second 
order into an orthogonal pair, consisting of vacuum and thermal noises which are 
orthogonal to the vacuum and connected by a reversible input-output (modular) 
filter. The time-continuous representation of the stationary linear filters requires a 
second order theory of quantum stochastic integration. 

The direct application of standard quantum stochastic integration with respect 
to the quantum white noise integrators is not possible for equilibrium (KMS) states 
because of the non-Markovian character of the corresponding standard quantum 
processes x and x. 

But it is possible to modify this approach for the spectral domain where the 
stationary quantum processes are (5-correlated. Using this approach in the second 
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order, we proved that they can be canonically decomposed into the superpositions 
of the vacuum pair of direct and reversed noises and can be obtained by a linear 
filtering from the standard vacuum noises x and x. 

In general the spectral quantum stochastic integrates need not to be adapted and 
homogeneous with respect to the frequency shifts and the treatment of quantum 
nonlinear filters and other transformations of quantum equilibrium noise requires 
non stationary and non adapted theory of quantum stochastic integration 
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